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The presented research was focused on a comparison between different means of obtaining a Nusselt number distribu- 
tion, in a situation where neither temperature nor heat flux density is constant. Two fundamentally different measure- 
ment techniques have been utilized, alongside a CFD simulation, in order to designate temperature distributions in a 
horizontal rod. Dry air under normal pressure, regarded as a perfect gas, was chosen as the working fluid, whereas 
the rod’s cross-section was restricted to a ring. In this scenario heat exchange between the rod and the fluid is 
driven predominantly by natural convection, with a slight impact of thermal radiation, particularly at temperatures 
approaching the top end of the available range. Temperature margins achieved at the heated end of the rod ranged 
from 60K up to 150K, resulting in local Rayleigh numbers falling in-between 6.0x10° and 2.6x10*. Reconstruc- 
tion of Nusselt numbers from a discrete temperature distribution was possible thanks to a dedicated method im- 
plemented using a Scilab script. A segregated, steady-state solver based on the SIMPLE scheme was utilized for 
the purpose of numerical simulations on the fluid side, whereas a heat conduction equation was solved over solid 
domain in the considered conjugated heat transfer problem. A corresponding set of empirical data has been ob- 
tained, using both resistance temperature detectors and a thermal imaging camera, both for the sake of numerical 


model validation and comparison of individual methods. The Nusselt numbers resulting from each approach were 


compared against values computed using available correlations valid for horizontal configuration. 
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Introduction 


Even though heat transfer from a straight rod is a ra- 
ther straightforward phenomenon, it bears a significant 
practical importance. This feature owes greatly to the 
similarities it shares with heat transfer from multiple 
kinds of fins, commonly found in thermal engineering 
practice. Not only the physical mechanisms are analog- 
ous in both cases, but also their mathematical descrip- 
tions are almost identical, at least at the most fundamen- 
tal level. Therefore, it appears to be solid testing grounds 
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for a number of techniques aimed at heat transfer en- 
hancement through surface area extension. 

Natural convection around rods and cylinders was in- 
vestigated by a number of researchers. The most popular 
studies concern temperature distribution and Nusselt 
number analysis. Atayilmaz [1] numerically and experi- 
mentally studied surface temperature distribution of a 
cylinder as a function of incination angle. Temperature 
distribution along a rod can be analyzed in various orien- 
tations as well: vertical (e.g. Basit et al. [2], Nagendra et 
al. [3]), horizontal (e.g. Atmane et al. [4], Liao and Lin [5]) 


NOMENCLATURE 


h heat transfer coefficient O=T-Ta temperature margin 

P perimeter O=6/45 dimensionless temperature 
m coefficient é=x/L dimensionless length 

Ta ambient temperature a dimensionless coefficient 


or as previously mentioned — at differing inclination 
(e.g. Heo and Chung [6], Al-Urabi and Salaman [7], 
Chand and Vir [8]). 

Arrangement of surroundings was proven to play a 
significant role in establishing a thermal equilibrium in 
scenarios where heat transfer is driven predominantly by 
natural convection. Liao and Lin [5] analyzed this phe- 
nomenon using a horizontal cylinder located in a square 
enclosure. Different types of enclosures were also ac- 
counted for by Wang et al. [9] and Lee et al. [10]. Most 
of the authors also took into consideration distribution of 
the Nusselt number. In the presented article heat transfer 
involving a horizontal rod is discussed, thus the 
semi-empirical formulas for Nusselt number developed 
by Morgan [11], Fand et al. [12] and Churchill and Chu 
[13] were chosen for comparison. 


Theory 


Analytic framework 


The basic theory of heat transfer in a straight, hori- 
zontal rod is based on the local equilibrium hypothesis 
followed by the energy conservation principle in a diffe- 
rential form. Consider a said rod with a constant heat flux 
delivered to one of its ends. Given that the rod's diameter 
remains small in comparison with its length, temperature 
should be approximately uniform over its cross-sections, 
rendering the case a 1-D problem. 

At any given distance x measured from the beginning 
of the rod there is an infinitesimal slice of length dx per- 
pendicular to the rod's axis. It experiences three distinc- 
tive heat fluxes on its boundaries: conductive flux at one 
end Q,, its counterpart at the other Qx+ax and a total out- 
going heat flux at its outer surface Qout. The energy con- 
servation principle states that: 


(Ò, = Ode) = dO oun =0 (1) 
Application of Fourier’s and Newton’s laws: 
Ò, = “kA, dO yy, = hPdxO (2) 
bs 
leads to the following form of the governing ordinary 
differential equation: 


2 
d°o (hP 
~ ma a 


which may be written in a dimensionless form, using 


č=x/L, O©=0/0, and m=hP/k,A, as: 
2 
Ze —(mLY @=0 (4) 
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The above equation possesses an analytic solution of 
the following form: 

Olé) = Cet + Ce" (5) 

In particular, a Dirichlet type boundary condition ap- 
plied at the heated end: ©ọ=1 and a Neumann type 
boundary condition at the other end: g,=0 lead to a rela- 
tively simple form of temperature distribution: 

cosh(mL(1 — &)) 
X= r a 
cosh(mL) 

The main deficiency of the above framework lays in 
omission of the radiative heat transfer. This problem may 
be overcome by perceiving the total heat transfer coeffi- 
cient h as a sum of convective and radiative contribu- 
tions: 


(6) 


h=hcthp (7) 
where: 
hp =e(T)o(T +T,)\T? +77) (8) 

Although this approach allows for retainment of the 
solution’s form of (5), it only does so under an assump- 
tion of a constant Ag. Keeping in mind that it constitutes a 
3" order polynomial of local temperature, this approach 
is only valid in case of a very small temperature differ- 
ence @ or unique materials characterized by an emissivity 
é(T) counteracting the temperature dependence. 

On the other hand, utilization of the total heat transfer 
coefficient defined in (7) is suitable both for application 
in a simple 1-D numerical scheme and elimination of 
radiative contribution from empirical results. Hence, nu- 
merical methods have been rendered a useful tool in 
analysis of the presented conjugate heat transfer problem. 
Furthermore, they allow to conduct an analysis in case of 
a more complex geometry, where experimental methods 
fail either due to limited resolution or accessibility. 


Average Nusselt number at varying T and q 


In order to obtain a Nusselt number distribution along 
the rod’s length from a discrete set of empirical data, a 
certain numerical method has been applied. It’s based on 
the 2" order central differencing scheme, although higher 
order numerics could be employed, given a denser dis- 


tribution of temperature. Using the previously introduced 
notation one gets: 


k k 
ae po = P=Nu-L6? (9) 
kA ky k; D k, 

where ¢° constitutes a coefficient describing the geo- 
metry of the rod and equals for a full rod: CAL yD’, 
whereas for a cylindrical pipe one obtains: 
CAE sel a), On the other hand, the dimension- 

less coefficient (mL) equals to the governing equation: 
ae ee 

(ney -4 

dx T-T A 
Equating the right-hand sides of (9) and (10) yields: 


(10) 


o (11) 


Let us assume, for the sole sake of simplicity, that a 
temperature distribution obtained from a set of linearly 
distributed values is given, yielding the spatial interval 
between each pair of neighbouring measurement points 
equal to Ax=L/n. Utilization of the 2nd order central dif- 
ferencing scheme implies that the second derivative at 
any point, excluding the boundaries, equals: 


AT Ty - 27, +T 


(12) 
dx? (Ax) 
which substituted into (11) leads to: 
ko n T,.,-27, +7, 
Nu; = S . i+] L i-l (13) 
ksi el T; -T1 
or using the dimensionless temperature ©: 
2 
Nu, - (Suits) (14) 
kif ©; 


At the boundaries the 2”? order right- or left-sided dif- 
ferencing scheme may be used instead. Finally, one ob- 
tains an average Nusselt number from: 


1 n 
Nu sy = f Nu(é)E = 1Y vu (15) 
0 i=l 


Experimental set-up 


All of the experiments have been conducted in a dar- 
kened, enclosed chamber, in order to efface the radiative 
impact of surrounding equipment and maintain steady 
conditions for natural convection. Dry air under atmos- 
pheric pressure was utilized as the working fluid. 


Experimental apparatus 


Experimental apparatus consists of a straight rod, a 
cartridge heater, a power supply system, a data acquisi- 
tion system, temperature sensors and appropriate har- 
nessing. Reduction of rod's thermal capacity leads to a 


faster achievement of a steady-state, thus an aluminium 
pipe has been used instead of a full rod. The pipe meas- 
ures 18mm in outer diameter, 16mm in inner diameter 
and 280mm in length. As shown in Figure 1, the neces- 
sary heat flux is supplied to the rod from the cartridge 
heater. Temperature measurements were performed using 
Pt100 resistance detectors glued to the outer surface of 
the rod and a single thermocouple located inside the hea- 
ter. The acquisition system collects the data with a fre- 
quency equal to 0.2Hz. The presented set-up was ac- 
companied by a thermal imaging camera, which provided 
an alternative measurement method. 


Temperature sensors 
Ti T, Ts T, T; Te T; Ts 


Heater 


DC Power 
Supply 


Fig.1 Scheme of the experimental apparatus 


Measurement’s accuracy 


The accuracy of the temperature measurements is de- 
termined mainly by the accuracy of the Pt100 sensors, 
data acquisition modules and sensors’ installation method. 
The accuracy of the Pt100 sensors depends on local tem- 
perature and was specified by their producer — under 
achieved conditions it does not exceed +0.4K. The data 
acquisition system is characterized by an accuracy of 
+0.15K. In order to account for the temperature gradient 
within the thermal baffle of glue which fixes sensors to 
the rod, an additional margin equal to -0.3K has been 
added. The sensors are distributed along the rod’s length in 
a linear fashion with positioning accuracy equal to +1 mm. 


Numerical set-up 


The numerical methods utilized in the presented re- 
search may be divided into two groups — a CFD simulation 
of conjugate heat transfer involving the rod and surround- 
ing fluid and an in-house utility based on equation (14), 
serving the purpose of post-processing temperature distri- 
butions obtained using individual means. An open-source 
solver used for the former part was chtMultiRegionSim- 
pleFoam, which constitutes a part of the OpenFOAM 3.0.1 
toolbox. 


Topology, mesh and boundary conditions 
In order to reduce the computational effort required to 


obtain meaningful results, a relatively simple topology 
has been chosen. It consists of a cuboid cut throughout 
by a cylinder representing the insides of the pipe, which 
were not taken into consideration in this analysis. The 
division into fluid and solid domains was performed on a 
supplied mesh using the topoSet utility, since individual 
regions are not singled out prior to a simulation run. 

The structured mesh generated using blockMesh utility 
consist of 24 blocks and 340032 hexaedral elements, 4032 
of which form the solid domain. A total of 30512 external 
faces compose 7 main boundaries — six walls of the cuboid 
and the pipe’s inner wall. All of them but the one at rod’s 
insulated end were assigned wall type BCs with no-slip for 
velocity field and either fixed temperature or zero normal 
gradient (rod’s inner wall) for temperature field. The last 
one was treated as a symmetry plane regarding all the va- 
riables. The rod was assigned boundary conditions analog- 
ous to those used to obtain analytic solution given by (6) — 
fixed temperature at the heated end and a symmetry type 
BC at the other end. 


Numerical schemes 


The chtMultiRegionSimpleFoam is a solver which 
combines the functionality of buoyantSimpleFoam and 
heatConductionFoam for the fluid and solid regions, re- 
spectively. BuoyantSimpleFoam, on the other hand, con- 
stitutes a SIMPLE based, steady-state solver for buoyant 
flows of compressible fluids. 

In order to take the rod’s irradiation into account, the 
viewFactor model was utilized with an appropriate 
coupling at the solid-fluid interface. Even though the 
convectively induced flow around the rod is purely la- 


minar, the k-s turbulence model was incorporated as well. 


It owes to the fact, that height of the domain is sufficient 
for corresponding Reynolds number to fall into the tran- 
sition regime. 

Among the most notable numerical schemes one 
should point out a bounded Gauss upwind scheme for 
most of the divergence terms and a Gauss linear scheme 
for the gradient terms. 


Data processing 


Nusselt number distributions and resulting average 
values were computed in all cases using a dedicated Sci- 
lab script. Transport properties of air were obtained sep- 
arately at each given point on the basis of a representa- 
tive temperature, constituting an arithmetic mean of local 
and ambient temperatures. The polynomial interpolation 
utilized for processing of the thermograms was based on 
the least squares method. 


Results and Discussion 


Three values of wattage were taken into consideration: 


17W, 33W and 50W. Approximately only one-third of the 
dissipated power is transported into the rod, due to losses 
at insulation connecting the rod to the frame. For this 
reason a Dirichlet type BC was supplied to the numerical 
calculations. 


Empirical results 


Two sets of data were gathered at each of the 3 power 
levels — one obtained using RTDs and the other resulting 
from thermal imaging. Even though both are discrete, 
there is huge gap between spatial resolution of these me- 
thods, ranging up to two orders of magnitude. There are 
only 8 sensors mounted to the rod as opposed to about 
300 pixels in the camera’s measurement area. Thus, the 
results originating from the optical method are presented 
using continuous lines. 

The temperature distributions obtained on empirical 
basis are presented in Figure 2. As anticipated, values 
from RTDs maintain a proportional difference and a good 
agreement with data from thermal imaging (TD, 
post-processed with an assumption of constant emissivity 
é=0.96. The data set characterized by e=0.60 is presented 
solely for comparison, since the value was used in a se- 
ries of numerical simulations. What’s troubling here is 
the systematic error at the 2"’ sensor from the heated end. 
It may be contributed to irradiation from the rest of the 
experimental apparatus, which also reaches a considera- 
bly high temperature. Nevertheless, conducting thrust- 
worthy measurements at this point was impossible, hence 
causing serious issues in context of computing an aver- 
age Nusselt number from this set of data. 
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Fig. 2 Temperature distributions along the rod’s length re- 
sulting from experiment 


Even though temperature distributions resulting from 
thermal imaging seem smooth, the spatial resolution of 
this method is relatively high, yielding differences be- 


tween subsequent values of the same order of magnitude 
as measurement’s accuracy. What follows is that the 2" 
derivative with respect to spatial coordinate, computed 
accordingly to (12), is highly unstable. In order to over- 
come this problem a 5" order polynomial interpolation 
has been utilized to compute the derivative and so Nu. 
Results of this approach, together with Nusselt numbers 
computed from the RTD measurements on the basis of 
(14), are presented in Figure 3. 

It is clearly visible that the utilized polynomials were 
not capable of an accurate designation of Nusselt number 
distribution, particularly in the cases of medium and high 
wattage. Although higher order polynomials brought sa- 
tisfactory results in terms of qualitative agreement over 
middle parts of individual distributions, they yielded very 
high discrepancies at the ends. Thus, they have been 
rendered inappropriate for designation of average Nusselt 
number, unless an arbitrary cut-off of distribution ends 
would be conducted. 

In spite of the described difficulties, approximation 
procedure returned satisfactory results in the 17W case, 
both in qualitative and quantitative terms. The results 
achieved under an assumption of emissivity «=0.96 
closely resemble those obtained using RTDs. This feature 
owes to a smoother initial temperature distribution ob- 
tained in this case, thus indicating that the problems with 
polynomials originate from the technical aspects of the 
experimental set-up, rather than from the method itself. 
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Fig. 3 Nusselt number distributions along the rod’s length 
resulting from experiment 


Numerical results 


The choice of boundary conditions for the numerical 
simulations reflects the conducted experimental work. 
Two different sets of data were gathered, with and with- 
out inclusion of a model responsible for radiative heat 
transfer. A flow pattern typical for the utilized configura- 
tion is presented in Figure 4, alongside rod’s temperature. 


The temperature distributions obtained under absence 
of radiation (NR) are milder than their counterparts re- 
sulting from simulations conducted using the viewFactor 
model (VF), as shown in Figure 5. It is, without a doubt, 
a consequence of the associated decrease of heat flux 
density, since the impact of radiation is usually consider- 
able at this range of temperature margin 0. 

Omission of radiative heat transfer resulted in an 
overestimation of the achieved temperature distribution. 
Regrettably, at emissivity <=0.96 temperature is incon- 
sistent with the RTD measurements, likewise. Artificial 
reduction of ¢ finally brought satisfactory agreement and 
an accurate estimation of average Nu. This feature possi- 
bly owes to some minor flaws in the numerical set-up or 
the viewFactor model itself. On the contrary, Nu distribu- 
tions obtained using this approach are relatively accurate. 
The characteristic teeth shown in Figure 6 result from 
agglomeration procedure utilized by the viewFactor 
model. 
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Fig. 4 Velocity streamlines originating from the mid-plane 
and rod’s temperature in the 50W NR case 


Thus, their amplitude and simultaneously length may 
be reduced at the expense of required computational effort. 


Comparison 


Four different types of semi-empirical correlations, 
valid for natural convection at Ra not exceeding 10°, 
were taken into consideration. The one developed by 
Morgan [11] implies that Nu is a function of Ra exclu- 
sively: 

Nu =0.48Ra°*> (16) 


Fand [12] et al. postulate that Pr should be incorpo- 
rated with a different exponent than Gr: 


Nu = 0.474Ra?” Pr?" (17) 
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Fig. 5 Temperature distributions along the rod’s length re- 
sulting from numerical simulations 
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Fig. 6 Nusselt number distributions along the rod’s length 
resulting from numerical simulations 


Oosthuizen assumes that Nu depends only on Gr, neg- 
lecting the impact of Pr completely: 
Nu =0.42Gr°” (18) 
Finally, a general correlation given by Churchill and Chu 
[13] was chosen for the comparative analysis: 


16 0.1685 


aA 
0.559 \16 
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Nu?’ =0.54+0.39| Ra 14 (19) 


Average Nusselt numbers computed on the basis of each 
of the above formulas, using temperature distributions ob- 
tained through the RTD measurements to find transport 
properties, are summarized in table 1, beside the empirical 
and numerical results. 

First of all, the average Nu resulting from thermal im- 
aging under an assumption of ¢=0.60 stand out as impre- 
cise, which is natural given their comparative purpose. 
Secondly, the other set of data obtained using the optical 
method is highly accurate in comparison with semi- 
empirical correlations, although it does not follow the 
expected pattern of Nu increase accordingly to the dissi- 
pated power. This could be a result of difficulties en- 
countered at the higher levels of wattage, associated with 
increasing roughness of temperature distribution. 

The values obtained using RTDs are defective, due to 
the problems with appropriate radiative insulation of the 
sensors. Nevertheless, they do not bring a severe unde- 
restimation. Hence, the presented method of Nu distribu- 
tion reproduction seems promising in context of small, 
discrete data sets. Furthermore, the results are expected 
to be improved, given an increase in spatial resolution. 

Finally, the numerical simulations constitute a precise 
and reliable source of knowledge about Nu distributions 
in the considered case. Not only they bring relatively 
stable results but also allow for a direct assessment of the 
radiative contribution to the overall heat transfer coeffi- 
cient. Moreover, they are suitable for complex geome- 
tries, where neither of the presented experimental tech- 
niques is applicable. On the other hand, they require both 
an initial validation and precise knowledge of emissivity 
of the considered surfaces as a function of temperature. 


Conclusions 


Even though direct sensor measurements constitute the 
most dependable method of obtaining a temperature dis- 
tribution in terms of absolute values, they suffer greatly 
from limited spatial resolution. Thermal imaging tech- 


Table 1 Comparison between empirical, numerical and semi-empirical results 


Nus empirical 


Nusy numerical 


Nusy semi-empirical 


Wattage 4 
F 
W) RTD  TIs=0.60  TIs=0.96 NR VF =0.60 VFe=0.96 Morgan 7 Oost-huizen — Chur-chill 
etal. 
16.6 3.106 2.044 5.085 3.830 5371 6.350 4.912 4.776 4.666 4.376 
332 3.319 2.346 4.715 4.097 5.762 6.771 5.443 5.293 5.174 4.851 
49.7 3.945 2.406 4377 4.200 6.003 7.064 5.652 5.495 5.374 5.039 


nique provides means to overcome this difficulty at the 
expense of achieving measurement accuracy of the same 
order of magnitude as difference between subsequent 
values, leading to an unstable 2" partial derivative with 
respect to coordinates. Fortunately, this issue could be 
mitigated using a polynomial approximation of the tem- 
perature distribution, which leads to a continuous and 
more accurate distribution of Nu. 

Numerical methods share some strengths and weak- 
nesses with thermal imaging, despite belonging to a com- 
pletely different group of techniques. The spatial accuracy 
of virtual measurements may be easily adjusted up to val- 
ues typical for the optical method. Furthermore, both me- 
thods require a priori knowledge of emissivity of the con- 
sidered surface as a function of temperature. Nevertheless, 
the obtained temperature distribution is rather smooth, al- 
lowing for a direct computation of Nu, instead of using a 
fitted curve. However, the numerical methods should not be 
applied without an appropriate validation, including a pa- 
rameter sensitivity analysis. 

All in all, thermal imaging technique maintains a com- 
promise between provided accuracy and time effort re- 
quired to set up. Under circumstances where its application 
becomes tedious, numerical methods may prove very handy. 
It should be kept in mind though, that they have to be em- 
ployed in a careful manner in order to bring quantitatively 
meaningful results. These methods seem superior to sensor 
measurements in context of an indirect identification of 
local Nusselt number. The latter still constitutes a sound 
alternative, particularly whenever thermal imaging is im- 
possible to conduct and numerical methods become unreli- 
able. 
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